Light-ion production from O, Si, Fe and Bi induced by 175 MeV quasi-monoenergetic 

neutrons 
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We have measured double-differential cross sections in the interaction of 175 MeV quasi- 
monoenergetic neutrons with O, Si, Fe and Bi. We have compared these results with model calcula- 
tions with INCL4.5-Abla07, MCNP6 and TALYS-1.2. We have also compared our data with PHITS 
calculations, where the pre-equilibrium stage of the reaction was accounted respectively using the 
JENDL/HE-2007 evaluated data library, the quantum molecular dynamics model (QMD) and a 
modified version of QMD (MQMD) to include a surface coalescence model. The most crucial aspect 
is the formation and emission of composite particles in the pre-equilibrium stage. 



I. INTRODUCTION 

Accelerator-driven incineration of nuclear waste may 
play an important role in future scenarios both in coun- 
tries investing in the nuclear renaissance and in countries 
deciding to phase out nuclear energy. However, develop- 
ment of accelerator-driven systems requires new and more 
reliable neutron-induced cross-section data in the inter- 
mediate region from 20 to 200 MeV. In this energy region 
the entire structural properties of the target nucleus may 
participate in the nuclear reaction, and the nuclear data 
community expressed the need to obtain benchmark ex- 
perimental data to improve present models and evaluated 
nuclear data libraries [l[. 

A series of neutron induced cross-section measurements 
at 96 MeV have been conducted at the The Svedberg Lab- 
oratory (TSL), Uppsala 0. However, experimental data 
are still scarce above 100 MeV. In this perspective, we 
have measured cross-sections for light-ion production in 
the interaction of 175 MeV quasi-monoenergetic neutrons 
(QMN) with C, O, Si, Fe, Bi and U. 

In this work, we compare O, Si, Fe and Bi experimen- 
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tal data with INCL4.5-Abla07 d Q, MCNP6 @ using 
CEM03.03 H, TALYS-1.2 and PHITS @ model cal- 
culations. While most models are able to reproduce the 
production of protons, the most crucial aspect is to de- 
scribe the formation and emission of composite particles 
in the pre-equilibrium stage. 



II. MATERIALS AND METHODS 

The experimental data presented in this work were 
measured with the Medley setup @ at TSL. Here a quasi- 
monoenergetic pulsed neutron beam is available, with 
maximum peak energy of 175 MeV. Time-of- flight infor- 
mation was used to reduce the contribution of the low- 
energy component of the neutron spectrum, as described 
in Ref. @. Model calculations presented in this paper 
were folded with the experimental neutron spectrum, for 
comparison with measured data. 

The Medley spectrometer is an array of eight three- 
elements AE-AE-E telescope detectors designed to mea- 
sure protons, deuterons, tritons, 3 He and ev-particles at 
eight angles in the laboratory system from 20° to 160°. 
The present configuration of Medley allows good parti- 
cle separation and at the same time a low identification 
threshold (2 MeV for protons) and a wide dynamic range 
(up to 170 MeV). Experimental details and data analysis 
procedures are described in Ref. @. 
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FIG. 1. Experimental deuteron production from Fe at 40° in 
comparison with model calculations. 
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FIG. 2. Experimental a particle production from Fe at 20° in 
comparison with model calculations 



Determination of the experimental cross sections was 
achieved measuring elastic (n,p) scattering from a CH2 
target at 20°; counts from other reaction targets where 
then normalized to the same neutron flux and the same 
number of scattering centers. 

III. RESULTS AND DISCUSSION 

To describe the pre-equilibrium emission of compos- 
ite light-ions, TALYS complements the two-component 
exciton model with direct-like mechanisms, according to 
the Kalbach systematics We observe how TALYS 

overestimates the production of deuterons (Fig. 1) and 
a particles (Fig. 2,3). We observe the same trend also in 
the pre-equilibrium production of other composite light- 
ions, at all emission angles. We may improve the predic- 
tive capabilities of TALYS by reducing the contribution 
of the nucleon transfer mechanism in the pre-equilibrium 
process (TALYS-modified) [ll| . This result was achieved 
introducing an energy dependent scaling on the cstrip pa- 
rameter in TALYS. 

However, since the TALYS code does not include mul- 
tiple pre-equilibrium emission of composite particles, the 
TALYS-modified model calculations underestimate the 
cross section at low energies in the pre-equilibrium stage 
of the reaction (Fig. 1, 2, 3). 

The INCL4.5-Abla07 calculations integrates the intra- 
nuclear cascade (INC) model 0] with the ABLA dc- 
excitation model [j]. For Fe and Bi targets, INCL-Abla07 
well describes the pre-equilibrium deuteron production 
(e.g. Fig. 1) and it is in good agreement with production 
of other composite light-ions (e.g. Fig. 2, 3). The INC 
model is generally used to describe reactions involving 
neutron energies above 200 MeV. However, our results 
show how the INC model may be applied to describes 
reactions with neutrons of lower energies. 

MCNP6 @ uses the Cascade-Exciton Model (CEM) [| 
as implemented in the code CEM03.03 as its default 
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FIG. 3. Experimental a particle production from Bi at 20° in 
comparison with model calculations 



event-generator. CEM assumes that nuclear reactions 
occur generally in three stages. The first stage is the 
INC. When the cascade stage of a reaction is complete, 
CEM uses the coalescence model to create high-energy d, 
t, 3 He, and a particles via final-state interactions among 
emitted cascade nucleons, already outside of the target. 
The second stage is an improved version of the modified 
exciton model of pre-equilibrium decay followed by the 
third equilibrium evaporation/fission stage of the reac- 
tion. But if the residual nuclei after the INC have atomic 
numbers with A<13, CEM uses the Fermi breakup model 
to calculate their further disintegration instead of using 
the pre-equilibrium and evaporation models. We see that 
MCNP6 calculations reproduce the emission of a parti- 
cles (Fig. 2, 3), but underestimates the forward emission 
of other light ions (e.g. Fig. I). 

Finally, we calculated double differential cross sections 
with the PHITS code Q, applying three different meth- 
ods to account for the pre-equilibrium stage of the reac- 
tions: using evaluated data from the JENDL/HE-2007 
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FIG. 4. Experimental deuteron production from O at 40° in 
comparison with model calculations 
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FIG. 5. Experimental a particle production from Si at 20° in 
comparison with model calculations 



library (T^ |. the quantum molecular d yna mics (QMD) 
model and a modified QMD (MQMD) [12| to include a 
surface coalescence model. In all three cases, the sta- 
tistical decay was calculated with the generalized evap- 
oration model. Whereas production of protons is well 
described by all calculations, the QMD model cannot 



account for the pre-equilibrium production of compos- 
ite ions (Fig. 4,5). The inclusion of a coalescence model 
in the QMD model improves the results for all compos- 
ite ions (e.g. Fig. 4, 5), however the underestimation of 
pre-equilibrium deuterons at forward angles (Fig. 1,4) 
seems to indicate the need to introduce direct-like reac- 
tion mechanisms in the pre-equilibrium stage. 

The discrepancies between PHITS calculations with 
data from the JENDL/HE-2007 library and experimen- 
tal production of deuterons (Fig. 4), tritons and 3 He, 
indicate the need to improve the evaluation of these 
reaction data. 

IV. CONCLUSIONS 

We have measured double-differential cross sections for 
several elements in the interaction with 175 MeV quasi 
monoenergetic neutrons, and we have compared our re- 
sults with model calculations. 

TALYS consistently overestimates the pre-equilibrium 
emission of all composite light-ions from Fe and Bi. We 
showed that reducing the contribution of direct-like pro- 
duction mechanisms, according to the Kalbach system- 
atica, improves TALYS prediction capabilities. MCNP6 
describes better the emission of a particles from Fe and 
Bi than of other composite ions, which are consistently 
underestimated in particular at forward emission angles. 
We have also seen that introducing a coalescence mecha- 
nism in the QMD model improves the predictive power of 
the PHITS code, and gives good agreement with a par- 
ticles data. Both MQMD and MCNP6 results seem to 
suggest the need to introduce direct-like mechanisms to 
account for the forward emission of deuterons, tritons 
and 3 He. Finally, INCL4.5-Abla07 describes fairly well 
the production of deuterons from Fe and Bi, but discrep- 
ancies emerge for other composite light-ions. 
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